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We derive the second-order hydrodynamic equation for reactive multi-component systems from the relativistic 
Boltzmann equation. In the reactive system, particles can change their species under the restriction of the 
imposed conservation laws during the collision process. Our derivation is based on the renormalization group 
(RG) method, in which the Boltzmann equation is solved in an organized perturbation method as faithfully as 
possible and possible secular terms are resummed away by a suitable setting of the initial value of the distribution 
function. The microscopic formulae of the relaxation times and the lengths are explicitly given as well as those 
of the transport coefficients for the reactive multi-component system. The resultant hydrodynamic equation with 
these formulae has nice properties that it satisfies the positivity of the entropy production rate and the Onsager’s 
reciprocal theorem, which ensure the validity of our derivation. 
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I. INTRODUCTION 

Under the highly extreme condition realized in the ultra- 
relativistic heavy ion collision experiments at the Relativistic 
Heavy Ion Collider (RHIC) at the Brookhaven National Labo¬ 
ratory (BNL) and the Large Hadron Collider (LHC) at CERN, 
the quark-gluon plasma (QGP) has been most probably cre¬ 
ated iSIl. To reveal the properties of QGP created there, 
the hydrodynamic model is utilized to study the time devel¬ 
opment of the extremely hot and dense matter. In particular, 
the large hadronic elliptic flow (^ 2 ) observed in the experi¬ 
ments indicates that QGP is described by the relativistic hy¬ 
drodynamic equation with a tiny viscosity, which may imply 
the creation of the strongly coupled matter iIs Ut^ : however, 
see also Refs. ids. Recent studies take into account such 
small viscous effects on the hydrodynamic expansion of the 
hot and dense matter to reveal the properties of QGP quantita¬ 
tively, where a relativistic dissipative hydrodynamic equation 
is an indispensable tool. A naive relativistic extension of the 
Navie-Stokes equation, however, has fundamental problems 
such as ambiguity of flow velocity, existence of unphysical in¬ 
stabilities, and lack of causality, and we need to introduce the 
second-order hydrodynamic equation to avoid such problems. 
Moreover, in the expanding object created in the experiments 
there are high density core which shows a nearly perfect fluid¬ 
ity and dilute corona in the peripheral region. In the latter part, 
more microscopic dynamics should be incorporated to the hy¬ 
drodynamics since the viscous effect is too large to apply the 
naive viscous hydrodynamics |@]. The second-order hydrody¬ 
namics actually include the mesoscopic dynamics and useful 
to analyze the hyd rody namic behavior of the spatially inho¬ 
mogeneous matter . 

The way of formulation of the second-order hydrodynam¬ 
ics is, however, controversial and many kinds of equation are 
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proposed. For instance, the Israel-Stewart equation, which is a 
kind of the second-order hydrodynamic equation widely used 
for the analyses of time development of the relativistic heavy 
ion collision or QGP, does not have established validity be¬ 
cause ambiguous assumptions are imposed to derive it. In fact, 
it has been shown that the solutions are different from those 
of the relativistic Boltzmann equation quantitatively ll20l - l^ . 

Recently, to eliminate the ambiguity in the derivation of hy¬ 
drodynamics and perform systematic formulation, the renor¬ 
malization group (RG) method ll33l - l^ has been applied to 
derive the second-order hydrodynamic equation for the single- 
component system both in the non-relativistic and relativistic 
cases ll29l - [^ as an extension of the first-order case ( 4 ^ : In 
this method, the Boltzmann equation is faithfully solved to 
obtain the distribution function around the local equilibrium 
state without imposing any ansatz in a perturbation theory and 
the secular terms are resummed away into the would-be inte¬ 
gral constants, which constitute the slow variables, i.e., the hy- 
drodynamical variables. It is also understood that, in the con¬ 
text of the reduction theory of dynamical systems 1^14^ . the 
hydrodynamical variables constitute the natural coordinates 
of the invariant/attractive manifold in the functional space 
spanned by the distribution function. Then the renormaliza¬ 
tion group equation ll^l^ or envelope equation ll^l^ gives 
the evolution equation of the hydrodynamical variables, i.e., 
the hydrodynami c eq uation, after an averaging of the distribu¬ 
tion function |43L l4^ . It has been shown that 

the resultant second-order relativistic hydrodynamic equation 
is causal and stable the microscopic expressions 

for the transport coefficients coincide with those obtained by 
Chapman-Enskog method, and those of the relaxation times 
are different from any other previous ones but allow physi¬ 
cally natural interpretations. Their numerical values are in¬ 
deed different from those by the moment method ll^ . 

In reality, the QGP is not a single-component system but 
reactive multi-component system composed of quarks, anti¬ 
quarks and gluons. Then it is imperative to derive the multi- 
component relativistic hydrodynamic equation for a realistic 
description of the QGP. In this paper, we apply the RG method 
to the derivation of the second-order relativistic hydrodynamic 
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equation for the reactive multi-component systems as an ex¬ 
tension of Ref. for the single-component system. In the 
reactive system, the species of the particles can be changed 
through the collision process in a way that the imposed conser¬ 
vation laws are satisfied; We restrict the collision to be a 2 to 2 
process for simplicity. We prove that the positive dehniteness 
of the entropy production rate and the Onsager’s reciprocal 
relation are readily satisfied by using our microscopic expres¬ 
sions of the transport coefficients. Furthermore, we show the 
causality and stability of the resultant equation in Ref. j^ . 
All these confirmation strongly suggest our formulation based 
on the RG method. 


Historically, Prakash et al. Ei first derived the second- 
order hydrodynamics for relativistic gas mixture as an exten¬ 
sion of the Israel-Stewart equation and investigated the prop¬ 
erties of hot hadronic matter. Then, there have been some 
attempts to derive the hydrodynamics for a mixture pOl - l^ . 
Among them, the work by Monnai and Hirano ISCT is of 
great importance in our point of view: They improved the 
Israel-Stewart’s fourteen moment method by additionaly im¬ 
posing two physically natural conditions, i.e., the positive¬ 
definiteness of the entropy production rate and the Onsager’s 
reciprocal relation for the transport coefficients to determine 
the functional form for the distribution function together with 
a careful order counting with respect to the Knudsen number. 
Thus they arrived at a multi-component relativistic hydrody¬ 
namic equation, which happens to have the same form as that 
of ours to be given in the present paper. It may imply that the 
form of the multi-component relativistic hydrodynamic equa¬ 
tion has been uniquely determined provided that the equation 
should satisfy the positivity of the entropy production rate and 
the reciprocal relations, irrespective of the microscopic forms 
of the transport coefficients. 


This paper is organized as follows: In Sec. [Ill we briefly 
review the basics of the reactive multi-component Boltzmann 
equation. In Sec. |III] we summarize theoretical foundations 
of the RG method. In Sec. IIVI we derive the second-order 
hydrodynamic equation for multi-component systems by the 
RG method. In Sec. m we discuss the basic properties of 
the resultant equation and the microscopic expressions of the 
relaxation times as well as the transport coefficients; a focus 
is put on the positivity of the entropy production rate and the 
reciprocal properties of the mutual transport coefficients. Sec¬ 
tion [VT] is devoted to a summary and concluding remarks. In 
AppendixiAl we work out to derive the functional forms of the 
excited modes utilizing the faithful solution of the Boltzmann 
equation. In Appendix iBl as an example, we apply the RG 
method to solve the Van del Pol equation . In Appendix O 
the derivation of the relaxation equation is performed explic¬ 
itly and the microscopic expressions of all the transport coef- 
hcients are given. 


In this paper, we use the natural unit, i.e., ?i = c = fee = 1, 
and the Minkowski metric = diag(-|-l, —1, —1, —1). 


II. MULTI-COMPONENT BOLTZMANN EQUATION 


In this section, we make a brief review of the relativistic 
transport theory based on the Boltzmann equation for reactive 
multi-component systems if^ . We consider a fluid composed 
of N species with M conserved currents. Let denote the 
charge of the /c-th species associated with the A-th conserved 
current, where k runs from 1 to and A from 1 to M. For in¬ 
stance, the electromagnetic U(l) charge of an electron is writ¬ 
ten as = —1. The conserved currents and the energy- 

momentum tensors are written as 


N 


= / ^PkPkfk,pdx), 

k=l 

N 

= Y^ / <iPkPkPkfk,pdx), 

1 _T ^ 


fc=l ‘ 


( 1 ) 

( 2 ) 


where fk,pk denotes the distribution function of the k- 
th species whose four-momentum is pk = {\/nik ~^Pk^Pk) 
with Pf, being the spatial components of the four momentum 

M 

Pk- 

Near-equilibrium behavior of the one-body distribution 
function may be described by the following multi-component 
relativistic Boltzmann equation ifs^ [3 : 

N 

Pk^kifkA^) = ^Cki[f]pdx), (3) 

Z=1 


with the collision integral 


Cki[fA{x) 



^ (1 T ^kfk,pk){X ^lfl,pi)y^ij\kl 

~ /fc.pfc//,pi (1 + Oi/i,Pi)(l + fj,Pj)y^kl\ij]- (4) 

Here, we have suppressed the arguments x, and abbreviated 
an integration measure as dp = d^p/[(27r)^p°]. Wki\ij is 
the transition probability due to the microscopic two-particle 
interaction with the energy-momentum conservation 

M 

yy^]\kl oc 5‘^{pi+ Pj - Pk - pi) , (5) 

A=l 

and the symmetry properties 

0 = Wy|/ci — Wp\ik, (6) 

tv 

0 = X! / dpzdpj [Wy|fc/ - Wfeiiy]. (7) 

The second relation is due to the unitarity of the scattering pro¬ 
cess ll3 . but note that Wij\ki ^ Wfcipj owing to the existence 
of a threshold of each reaction, ak represents the quantum sta¬ 
tistical effect, i.e., ak = -fl for boson, ak = —1 for fermion. 
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FIG. 1. Schematic figure of a multi-component collision. Though 
only binary collision occur, kinds of particles may change during 
each scattering process with charges conserved, i.e., qf -\- qf = 
Qk + ■ 


and afe = 0 for the classical Boltzmann gas. In the colli¬ 
sion processes described by this collision integral, we consider 
only binary collisions which conserve the charge before and 
after collisions. We, however, allow reactions, i.e., the parti¬ 
cle numbers of each component do not need to be conserved 
in a collision process. The allowance of the change of the total 
number of particles during each scattering process is of phys¬ 
ical importance but it is left as a future project. When only 
elastic scattering are allowed, we further impose = Sk^A 
and yVfc/|y = (1 - Ski/2){6kzSij + SkjSii)Wki, where Wu is 
the transition probability for the elastic scattering. 

On account of the symmetry properties (|6]) and O, the fol- 
lowing identity for the collision integral for an arbitrary vector 
fk,pk 1155n can be shown, 


In the Boltzmann theory, the entropy current may be de¬ 
fined by 


N 


= -J2 ^p>^Pk 


k^l ' 


fk,pk In fk,p 


(1 + akfk,pk ) ln(l + akfk,pk) 


ak 


(13) 


The entropy-production rate reads 


k,l=l'' ^ 


, 1 T ^kfk,p 


, (14) 


because of Eq. ©. One sees that the entropy production rate 
vanishes only if ln(/fc_p^/(I + o,kfk,pk)) is a collision invari¬ 
ant, i.e., ln(/p/(l + a/p)) =^fc.p, = <lkO^A{x) + 

pk-j3^{x). One thus finds toIImI] that the entropy-conserving 
distribution function can be parametrized as 


fk, 


Pk 


exp 


(pk-u 

\ ^ 


_ ^eq 


M 

E 

^=1 




- ak 


(15) 


E / ‘^PkV’k ,Pk ^kl [f]pk 

k,l=l 4 

1 .k, r 

= ^ E / dpidpjdpfedpi 

i,j,k,i=i-^ 

^ i^k,pk 4 ” ^l^pi ^i^Pi ^j,Pj )'^ij\kl 
^ /z,Pi /j,Pj (1 4” akfk.pk ) (1 4“ a^//p^). (8) 

Substituting {qk,Pk) iti^o </5fc,ps, in Eq. ©, we find that 
(Qk’Pk) collision invariants satisfying 


tv 

E 7fc / dpfeC'fei[/]p, = 0, 

(9) 



TV 

E / ‘^PkPkCkl[f\p^ = 0, 

(10) 


which mean the charge and energy-momentum conservation 
in the collision process, respectively. These two equations 
leads to the balance equations 

dpN^ = 0, (11) 

= 0 . ( 12 ) 


which is identified with the local equilibrium distribution 
function by interpreting T, pA, and as the local tem¬ 
perature, chemical potentials associated with the conserved 
charges, and flow velocity, respectively, with the normaliza¬ 
tion 

= 1. (16) 

We see that the collision integral identically vanishes for the 

local equilibrium distribution as 

Cki[r\, = o, ( 17 ) 

owning to the detailed balance 

tv 

^ ^ / ^Pi^Pj[fi,Pifj,Pj{^ 4- akfk,pk){^ 4- aifi^pi)yVij^ki 
~ /fc,Pfc/i,Pi(l 4- ai/i,pi)(l 4- k’-jfpPj 

N 

= ^ ^ / ^Pi^PjyVij\ki[fi,pifj,pji^~^akfk,pk){^~^aifi^pi) 

~ fk,pkfi,pi (1 4- ai/i,pi)(l 4- k'‘jfj,Pj)] 

= 0, (18) 


It should be noted, however, that any dynamical properties are 
not contained in these equations unless the evolution of fk,pk 
has been obtained as a solution to Eq. ©. 


where the symmetry property © and the conservations of 
energy-momentum and charges © have been used in the first 
and second equality, respectively. 
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III. FOUNDATION OF RG METHOD 


In this section, we explain the renormalization group 
method used to identify slow variables of the original dy¬ 
namics and obtain the reduced dynamics in terms of the slow 
variables JHHHl]. Here, we explain the method following 
the discussion in ll^ [dd]] . Let us consider the following n- 
dimensional differential equation; 

^ = (19) 

We denote W (t) for the exact solution of Eq. ( fT^ though it 
is unspecified yet. By setting the exact solution W {t) as the 
initial condition, we solve Eq. ( fT9l l perturbatively to obtain 

dX{t,to^W{to)) ^ p^x^t;to,Wito)),t). (20) 

The basic concept of the renormalization group method is that 
the true solution of Eq. ( fT9b is independent of the initial time: 

X(f; <0, W{to)) = X{t; f', W{Q), (21) 

which, by taking the limit to, leads to the RG equation 


dX{t;to,Wito)) 
dto 

Since the perturbative solutions X{t;to, X{to)) and 
X{t-,t'o, X{t'o)) are only valid around t ^ to and t ^ fg, 
respectively, the natural condition for f is to < f < tg which 
gives the condition to — t after taking the limit tg ^ tg. 
Now, Eq. (l22l) may be interpreted as a condition to construct 
the envelope of the family of curves which is represented 
by X{ty^, X(tQ)) with tg being a parameter characterizing 
curves l3^ . By solving the RG equation, we can write 
down the RG improved solution, or envelope, as 

XG{t) = X{t;to,W{to))\to=t. (23) 

First, we check that the global solution (|2^ satisfies the 
original differential equation ( fT9] l. Substituting Eq. ( |2^ into 
the left-hand side of Eq. ( fT9] l. we have 



dX{t-,to = t,W{t)) 
dt 


dX{t;to,W{to)) 

dt 


to=t 


dXit-,to,W{to)) 


dfg 



Since X{t;to,W{to)) satisfies ( |20] | and because of the RG 
equation (ESl l. Eq. (l23t solves Eq. (fT^ : 




F{XG,t). 


(25) 


Next, we show the converse, namely, that X{t; to, W(to)) 
satisfies the RG equation (1221) . provided that Xq (t) given in 
Eq. (|2^ solves the original differential equation (fT^ : 


FiXG,t) 


dXG 

dt 


dX{t-,to,W{to)) 


dt 


to=t 


dX{t-,to,Wito)) 

dto 


to=t 

(26) 


Since X{t-,to,W{to)) satisfies Eq. (l20l l. the satisfaction of 
the RG equation (l22l l is readily confirmed. 

Therefore, the RG equation (|22| and the original differen¬ 
tial equation with the global solution inserted ( l20b are equiv¬ 
alent. In the next section, we solve the Boltzmann equation 
following the latter method although some elaboration is nec¬ 
essary in extending the invariant manifold so as to incorporat¬ 
ing the excited modes 1^ . In Appendix iBl the equivalence 
of these two methods are demonstrated by explicitly working 
out an example with a limit cycle, although an incorporation 
of the excited modes is not made. 


IV. DERIVATION OE SECOND-ORDER 
HYDRODYNAMICS EOR MULTI-COMPONENT SYSTEMS 

The derivation of the hydrodynamic equation by the RG 
method explained in the last section consists of the following 
three steps: 

1. Solve the relativistic Boltzmann equation perturbatively 
with respect to small spatial inhomogeneity. 

2. Improve the perturbative solution to obtain the global 
solution with the renormalization group. 

3. Substitute the global solution into the Boltzmann equa¬ 
tion and construct the moment equations, which give 
the hydrodynamic equation. 

In this section, we utilize the RG method to drive the 
second-order hydrodynamic equation for multi-component 
systems from the relativistic Boltzmann equation ( l27l i. Then, 
we discuss some properties of the resultant equation, focusing 
on the entropy production rate. We shall see that the resultant 
equation readily satisfies the positivity of the entropy produc¬ 
tion rate and the Onager’s reciprocal relation, which is the 
general properties of the non-equilibrium phenomena. 


A. Solving the Boltzmann equation to extract mesoscopic 
dynamics 

1. An organized perturbation theory 

To solve the Boltzmann equation perturbatively with re¬ 
spect to spatial inhomogeneity, which is assumed to be small 


dt 
















5 


since the system is in the hydrodynamic regime, we rewrite 
Eq. Q as 


d 1 

T^fk,pdr,cr) = - ^Cki[f]pdT,cr) 

- ^——Pk-^h,pAT,(y)- (27) 

Pk-u oai^ 

where e is a bookkeeping parameter of the small spatial in¬ 
homogeneity corresponding to the Knudsen number = (mean 
free path) / (characteristic macroscopic length scale). After 
deriving the hydrodynamics, we shall put back e = 1. Here, 
we have introduced the temporal and spatial derivatives 

^ (28) 

respectively, is interpreted as the flow velocity satisfying 
UpU^ = 1, and — u^u''. The present choice of 

the coordinate system leads to the relativistic hydrodynamics 
in the Landau-Lifshitz frame or the energy frame ^56l] . From 
now on, we suppress the variable tr. We expand the distribu¬ 
tion function in the following perturbation series with respect 
to e: 

f{T;To) 

= + e/(^)(T;ro) -f eV^^^(T;ro) + O(e^). 

(29) 


Eq. ( |29] | is understood as the expansion around the equilib¬ 
rium distribution function and with i > 1 represents 
the non-equilibrium correction, which gives rise to dissipative 
effects. 

The hrst- and second-order equations read 


A 

A 


(35) 

+ rq/eqiT(T-ro), (36) 


with the definitions of the following quantities 


= -ifklJklJ-'^P'^ • V/il, 




k,pk\m,qm',ii,rn. 


, req req n -1 

k,pkJ k,pk> 


(37) 

(38) 

(39) 


,52 


Sfm ,qm^fn.rn \Pk * 


1 

— [f]k,pk 


^(1) ^(1) 
A.m,qTn Tn,rn 


/=/*■> 


(40) 


and the linearized collision operator 


Correspondingly, we expand the initial value as 
/(t = To]To) 

= /(to) 

= f^°\To) + ef^^^To) + e2/(2)(To) + 0{e^). (30) 

which is assumed to be an exact solution of Eq. (l27l l. Sub¬ 
stituting them into Eq. ( l27b . we obtain a series of equations, 
which we solve order by order. 

The zeroth order equation with respect to e reads 


Bk, 


Pk;m,qm 


— I fsq req n -1 

~ \J k,pkJ k,pk> 


X 


5 

^fm,qm 


\Pk • U 



feq 

lm,qm 


xeq 

Jm,qm 





1 

‘2pk ■ u 


dpldp^dpJyVij\kl 


Jeq peq peq 
i,piJ j,pj 
peq 


^ \^^km^{Pk Qm) d~ Slm^{Pl Qm) 

^im^iPi Qm) ^jm^iPj 9m)}; (41) 


AO) 

QpJ k,pk 


{t-,to) 


1 

Pk ■ U 


N 


'^Cki[f^%^ (r;To). 
1=1 


( 31 ) 


Since we are now considering a near-equilibrium system, 
the zeroth-order solution should satisfies 


^ f(o) 1 ^ n 

-^fk,p,ir]To) = 0, 


(32) 


With respect to the inner product defined as 

N 

(V'-x) = 

k=l 

for arbitrary vectors 'il}k,pk Xk,pk^ the linearized colli¬ 
sion operator Lk,p^-^m,qm found to be self-adjoint and semi¬ 
negative definite; 


dpkiPk 


u)fr 




k,pkJ k,pk 


(42) 


which leads to 


- E to) =0 (33) 

Pk ■ u 

Therefore, from Eq. (fTTl i the local equilibrium distribution 
function is found to be the zeroth-order solution: 


(34) 


{i’,Lx) = {Lip,x), ( 43 ) 

{ii, LiP) < 0. (44) 

The operator L has the M -f 4 eigenvectors belonging to the 
zero eigenvalue, i.e., 

a ^ [ Pfc’ a = A = 0,l,2,3, 

" I qt a = 4l + 3 = 4,.-. ,M + 3, 


f^^\T]To) = r{To). 


( 45 ) 
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which are the collision invariants and satisfy 

(46) 

We call pj, the zero modes. 

The Pp-space is defined as the space spanned by the zero 
modes of the linearized collision operator, which govern the 
slow dynamics of the solution, and the Qo-space is defined as 
the complemental space of the Pg-space in the solution space 
of the relativistic Boltzmann equation, which are spanned by 
excited modes. The projection operator onto the Pg-space and 
the Qo-space are given by 

k,p, = vlp, V~p , V’), (47) 

Qo = 1 - Po, (48) 

where is the inverse matrix of the the Po-space metric 
matrix 77 “^ defined by 

(49) 


We further divide the Qg-space into the Pi- and Qi- 
space: Pi-space is spanned by the doublet modes IBIT . 
which consists of the sets of ( 11 ^ ^^, JA,k,p,^Kl,) and 

, while the Qi- 

space is the complement to the sum of the Pq and Pi spaces; 
Here, If, J^, and are found to be the microscopic 
representations of the bulk viscosity, charge diffusion, and 
stress tensor, respectively, whose definitions are given in Ap¬ 
pendix. 13 The correspondent projection operators are de¬ 
noted by Pi and Qi, respectively. 

Equations (iTSl l and (l36l l can be solved much the same way 
as in El using a vector representation: We refer to El for 
the details. The first- and second-order solutions take the fol¬ 
lowing forms: 






(50) 


p\r) = nnUr - To)Po + (r - To)gis)-^ Pig{s)Qo 


- (1 +[t-TQ) d/ds)Qig{s)Qo K(s) 

with the initial values 


s =0 


, (51) 


/(i) = 

= -nnQig{s)QoK{s) 


(52) 

(53) 


After a lengthy calculation, which is worked out in App. |3 
we obtain the initial value T' for the first-order equation: 




[L-ifl] 


k,Pk 


n 


-f 


(n,L-in) 

M 

A.B^l 

L-I k,pk 


Jb,p 


where h = {e + P)/n is the enthalpy per particle and the 
matrix notation is introduced as: 

(>, L-^Jp)ab = {J^A^ L-\JbA. (55) 

and the (A, i?)-component of its inverse matrix is denoted by 
{JP, L~^ Jp{s))^g. The following would-be 3M -|- 6 integral 
constants have been introduced as the amplitudes of the nine 
vector fields; 


n(cr;ro), J^(cr;ro), 7r'"‘'(cr; tq). (56) 

We note that both and ttP'' are transverse and satisfy the 
equalities 

= AP''Ja,., (57) 

TtP‘' = AP'^P’^TTp^, (58) 

where AP'^p'^ is a traceless symmetric projection operator 

-f AP^^A'^P - 2/2,AP''AP'^). (59) 

We note that H, J^, and ■kP'^ will be interpreted as the bulk 
pressure, thermal flux, and stress tensor, respectively. We have 
also introduced a “propagator” defined by 

g{s) = {L - d/ds)~^, (60) 

which represents the temporal non-locality, i.e., non- 
Markovian effect. 

Summing up the perturbative solutions up to the second or¬ 
der with respect to e, we have the full expression of the initial 
value and the perturbative solution around t ~ rg to the sec¬ 
ond order: 


hr) = r + enn (1 + (r - t,)L)^ + (t - rgjFg 


_j_ g2 jeq jeq 


(r - Tg)Pg -f (t - Tg)C/(s) Pl(7(s)(3g 


- {1 + {t - To)d/ds)Qig{s)Qo K{s) 


s =0 


f = n + - e^nrQig{s)QoK{s) 


(61) 


s =0 

(62) 


which has completely the same structure as the one in the 
single-component case in the vector notation E3] . 


2 . RG improvement and moment equations 


The perturbative solution dSB possesses secular terms, 
which become divergent and clearly breaks down the pertur¬ 
bation as |t — Tg I gets larger. Improvement with the RG equa¬ 
tion allows us to avoid such secular divergences and obtain the 
global solution E3l - [35l4lll^l^ . By applying the RG equa¬ 
tion to the local solution dMT l. we obtain the rg-dependence of 
the would-be integral constants, which validate the approxi¬ 
mate solution in a global domain: 


1 fk,pk{r,o",TQ) 
drg 


= 0 . 


(54) 


To = T 


(63) 
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The RG equation (l6Tt indeed gives the equation of motion 
governing the dynamics of the would-be 4M + 10 integral 
constants T(cr; r), fj,a{cr; r), t), n(cr; r), J^(cr; r), and 

t). The global solution can be obtained as the initial 

value I 


leads to 


N 






{vk ■ u)— + epfc • V 


/ (t) = f (to = t) 

= n{i+er^) 

- e^nnQ,g{s)QoKis) 


s=0 


(64) 




= 0 + O{e^), (70) 

and the equation reduced by the excited modes reads 


To=T 


where the solution to Eq. (l6Tt is to be inserted. We have de¬ 
rived the mesoscopic dynamics of the relativistic Boltzmann 
equation (l27l l in the form of the pair of Eqs. (l6^ and (l64l l. Let 
us remark on the global solution obtained here; 


k—1 ^ 




r ^ d 

[Pk ■ u)— + epfe • V 


fG ^ yeq g2y(2) 


with 


f(2) ^-f<if-<iQ,g(s)QoKis) 


s=0 


(65) 


( 66 ) 


= e(L-i(n,J^,7r'^'^),Lvl/) 

+ e^^{L-\fl, Z,7T^n, B[^, 4/]) + 0(e3). 


(71) 


Without loss of generality, we can suppose that 'k contains 
no zero modes, because such zero modes can be eliminated 
by the redefinition of the zeroth-order initial value. It means 
that the possible existence of the zero modes in 4^ would be 
renormalized into the local temperature, chemical potential. 


Eqs. (iTOl i and (1711 1 give the equations of motion governing T, 
Pa, uk-, If, J^, and TTk" in the Pp-space and Pi-space, respec¬ 
tively. 


B. Hydrodynamic equation 

Now, we insert e = 1 in Eqs. (ITOl l and (fTTl i. Then, Eq. (ITOl l 


and flow velocity, where the local frame is defined and fixed 

is reduced to the following form 


by the flow velocity uk(a; tq). This is a kind of the matching 
condition. We note that the expression (|66]| tells us that 
does not include Pi- modes. Summarizing the above, we have 
the following orthogonality conditions 

with 

^M'^hydro ~ 6, 

(72) 



N 


0 = (+“,4-) = (+“, m 

hydro 

= Y ^PkPkPipkfZk^^ + fZk 

L, _1 ^ 

'^k,pk ) 



K —i 

J euku'' - {P + U)Ak'' +Trk\ 

a = V, 

= {L-\llJk^7Tk''),(f-^f-‘i)-^fi^)). (68) 


\ HAUk+ 

a = A + 3. 

Although, in the usual RG method, the RG equation (l63Tl 



(73) 


describe the slow dynamics of the would-be integral con¬ 
stants, the insertion of the global solution (l64l i into the 
Boltzmann equation (O give the equivalent dynamics for the 
would-be integral constants. Substituting Eq. ( l64l i into Eq. (O, 
we have 

= e{pk ■ 'k')fZi_fZk^k,Pk',m,qm^m,qm 

+ e\pk • 

+ e'^iPk-u)fllJ^l^^Bk,p,[f^^\f^%, (69) 

which is further reduced to simpler forms by taking the in¬ 
ner product with the zero modes and the excited modes 

[L“^(n, J^, ^ used in the definition of '^k,pk, re¬ 

spectively. The averaging inner product with the zero modes 


Here, n, e, and P denote the particle-number density, internal 
energy, and pressure, respectively, whose microscopic repre¬ 
sentations are given by 


UA = Ufj,N^ = Y^k dpkipk ■ u)fk,pk, (74) 

k=l •' 

N 

e = u^UuTk’' = E / ‘^PkiPk ■ u)^fk,pk, (75) 

k=l 

1 1 ^ f 

P = = - 3 E / ’ (76) 




We remark that Eq. (|72] t is nothing but the balance equa- 
tions and and can be identified with the energy- 

momentum tensor Tk'' and charge current given by Eqs. 
(|2li and ([T]i, respectively. 
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After a straightforward but lengthy manipulation worked 

I 


out in Appendix. O we can reduce Eq. (fTTI i into the following 
relaxation equations: 


M 


n = -C0-rn^n-^f“^V- 


M M 

.(1)^7. Vjprp^ Sr 


Ja,pVP^ + 


H- ^ ^ k\jj Ja.p^^T + ^ ^ rvYij . rji 


A=1 

M 


A.B^l 


+ ^tiniin^ + ^ ^ ^ujjJa’^b.p + ^iitt 

A,B=1 

M rp2 AI 

1 ^uf^B 


tj-TT '^pat 


d 


B^l B^l ^ 

M MM 

+ + J2 + J2 + E + Kfj^^JB,pU^^'^ 


K, 


(1)^ 

J-JT 

M 


S=1 

M 


B=1 


S=1 


irf^PVpT + J2 


S=1 

M 


Mab^pp 

JTT H rp 


+ E + E 


'7r‘ 

p'‘ 5 

M 


pp 


S=1 


B=1 


TTp, - 




a—1 


(77) 


(78) 


M M 


A^l 

M 


A,B^1 


+ 5^n.n7r'^"+ E b^fjJA^B 


(79) 


A,B=1 


where A pa is a traceless symmetric tensor. 

The scalar expansion 6 = V ■ u, the shear tensor = 
A^''P'^Vpa, the vorticity term = \{y^u'' — V^u'^), and 
many coefficients have been introduced, where explicit dehni- 
tions are given in App.lCl You also hnd the relaxation equation 
with e explicitly shown in App.lCl 

Here, we write down the resultant microscopic representa¬ 
tions of the transport coefficients C,, Xab, p, and relax¬ 
ation times Tn, , and r^, as follows: 

C = -^(n,L-in) 

ds(n(o),n(s)), (80) 

Xab = -^{JaA~^Jb,p) 

= -^l ds(J^(0),JB.^(s)), (81) 



where we have defined the “time-evolved” vectors by 


M 


-E 


J k,pk-,m,qm 


(83) 


Here, we note the symmetry property: Xab = Xba- The 
relaxation times have the expressions given by 


^ (n,£-^ri) _ Cd5g{n(o),n(5)) 

(n,T-iri) /“ds(n(o),n(s)) ’ 

M 

^ - E {L-^J^,L-^Jp)^A'^^,L-^Jp)cB 

c=i 
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M 


= E 

C=1 


dss(J^(0), J^(s)} 


AC 


ds(J^(0),J^(s)) 


(85) 


CB 


Tn = - 


{TTP'^,L-^npa) /g“ds(7r^‘^(0),7r^^(s)) 


( 86 ) 


We note that the cross correlations are generated by Xab, 
and so on for A ^ B. The existence of the cross corre¬ 
lations is characteristic feature of the multi-component fluid. 
Owing to the terms including such transport coefficients, a 
diffusion of some charge influences that of the others. 

It is noteworthy that C, Xaa, '^n, and can be 
proved to be positive definite as a result of the semi-negative 
definiteness of the linearized collision operator L as given in 
Eq. dSl). 

Here, we prove the positive dehniteness of the bulk vis¬ 
cosity C and the proof for the other coefficients are straight¬ 
forward. Since L~^ is symmetric and negative dehnite when 
acting on the excited modes, there exists a real (lower trian¬ 
gular) matrix U, such that L~^ — —U*U, which is called 
Cholesky decomposition. By using the Cholesky decomposi¬ 
tion, we have 

(n, L-^n) = -(n, u*uii) = -(c/n, ufi) < o (87) 

which show the positive definiteness of the bulk viscosity dSOl) . 

Let us clarify the difference between our result and the pre¬ 
vious attempts 14^1 - 1^ to derive the second-order hydrody¬ 
namic equation for relativistic multi-component systems. In 
dm, the Israel-Stewart theory is simply extended to multi- 
component systems and hence the resultant equation is not 
free from the drawbacks that the Israel-Stewart equation pos¬ 
sessed, even apart from the shortcomings that the reactive ef¬ 
fects are not included. In [1^] . though they also used the mo¬ 
ment method, they obtained new second-order terms by per¬ 
forming the careful order counting with respect to the Knud- 
sen number and indeed the resultant equation is consistent 
with ours. However, since the collision integral of the Boltz¬ 
mann equation was not specihed in their derivation of the hy¬ 
drodynamics, explicit expressions of the transport coefficients 
were not given. In addition, the Onsager’s reciprocal theorem 
was imposed a priori to determine the ansatz for the distri¬ 
bution function. In contrast, the microscopic forms of all the 
transport coefficients are given explicitly in terms of the col¬ 
lision operator in the present work. Furthermore, it is found 
that the very microscopic expressions of the transport coeffi¬ 
cients have so nice properties that they readily lead to not only 
the positivity of the entropy production rate but also the On¬ 
sager’s reciprocal theorem, which we shall prove in the next 
section. 


V. DISCUSSIONS 

We now examine the basic properties of the resultant hy¬ 
drodynamic equations (|72]) and (IC61l) - (IC631 l focusing on the 


entropy production rate. Since the one-body distribution func¬ 
tion (l64l l is an approximate solution of the Boltzmann equa¬ 
tion valid up to in a global domain, the positivity of the 
entropy production rate and the Onsager’s reciprocal theorem 
must be satisfied as the Boltzmann theory 1531] . In this sec¬ 
tion, we check that these properties are indeed satished by 
Eq. (|64] | up to e, since the terms of O(e^) do not contribute to 
the hydrodynamic equation due to the orthogonality relation 
dlTli and dSl. 


A. Positivity of entropy production rate 

The one-body distribution function ( l64b up to hrst-order of 
e is converted to the following form: 


fk,pk (T) — t) 


exp 


Pk • u 
T 


M 





Oie% 

( 88 ) 


where 'd>k,pk is given by Eq. (IA20I) . 

Inserting Eq. ( l88l l into Eq. (fT4l) . we have the entropy pro¬ 
duction rate in the following form. 


d„. 


k^l \ 


M 


-Pk^p E ^k^ - ^Pk9p'fk,pk 1 + O(e^) 


A=1 


T 


A=1 ' ' 


0{e^) 


I J_ 


‘ fiiy 


\TC 2Tr] 

,2 M \ 

E {^-^)ABJpB,A+0{e^), 


(89) 


A,B=1 


where the relaxation equations in the leading order is given by 


n = -C0-rn^n 

OT 


0(e), 


(90) 


M 


T2 




M 




+ 0{e), (91) 

= 277(7^" - + 0{e). (92) 

or 

In the second equality of Eq. ( l89b . we have used the following 













10 


transformation 

N 


ST'^'^ = E / ^PPkPUZA^'^^P'^ = + 7^'^^ 


k^l 


N 


(93) 


SK = ll^k j^PPtfllJll,^k,^u = Jl, (94) 

fc=i •' 


and 


where Ji is dissipative currents. Generally, the sum is taken 
over any tensor structure of Ji and Xi. In the linear approxi¬ 
mation, Ji is given by the following linear combination; 

J, = Y,JvX,, (98) 

3 

where all Xj in the sum have the same tensor structure as Ji 
and the coefficients are scalar quantities. The Onsager’s 
theorem asserts the symmetry property of the coefficients; 


= nw|:n + |5 E 


TQ dr 

'_L±. 

2Tr] dr 




A,B,C 
f 0(e). 


(95) 


In addition, the decomposition of the derivative, 0^ = 

Uf^d/dr + eV^, and the orthogonality between the fluid ve¬ 
locity and the dissipative contributions to the energy mo¬ 
mentum tensor and the particle current, and SN^, have 
been used. 

The viscosities ^ and rj are positive quantities as remarked 
in the last section. Positive definiteness of the third term of 
Eq. ( l89b is proved as follows; 


7»j=7i*- (99) 

From the expression of the entropy production rate ( |89] |- 
(l92b . one readily confirms that the relation ( l99b is satis¬ 
fied in the present case by identifying Ji = J'^, Xi = 
-(l/ft.)V^(pA/T’), and -fij = -XabT'^, because Xab = 
Xba from the microscopic expression of the transport co¬ 
efficients ( ISTT i. Since this relations are originated from the 
time-reversal symmetry of the underlying microscopic the¬ 
ory, this confirmation supports the validity of our treatment of 
the multi-component fluid; Note that the collision integral is 
given in terms of the microscopic transition rates that respect 
the time-reversal symmetry. 


VI. SUMMARY AND CONCLUDING REMARKS 


,2 M 

A,B = 1 

.2 M 

A.B=1 

.2 M 

= -^ E {JA,L-^jBAi^-^jnA{x-^Ji.)B 

A.B^l 

,2 I M M \ 

= ^ ( E UJ'^iX-^J^A, E UJbA>^~"J>^)b ) 

\yl=l B = 1 / 

> 0, (96) 

where we have substituted dSD) in the first equality and used 
the Cholesky decomposition in the second equality. As a re¬ 
sult, Eq. ( l89l l explicitly shows the positivity of the entropy 
production rate. 

B. Onsager’s reciprocal theorem 

Next, we consider the Onsager’s reciprocal theorem, which 
is originated from the time-reversal symmetry of the underly¬ 
ing microscopic theory 

The statement of the theorem is as follows. Near equilib¬ 
rium, the entropy-production rate takes the following form; 

^ J,X,. (97) 


In this paper, we have derived the second-order hydro- 
dynamic equation for the reactive multi-component systems 
from the relativistic Boltzmann equation based on the RG 
method. 

We have treated the system composed of N species with 
M conserved currents. In this case, the hydrodynamic vari¬ 
ables are composed of M -f 4 moments - local temperature, 
fluid velocity, and local chemical potentials conjugate to each 
conserved charges, which correspond to the conservation of 
energy, momentum, and M conserved currents. Similarly, the 
excited Pi-space that is incorporated in the second-order hy¬ 
drodynamics is spanned by 3M + 6 vectors because there 
are dissipative currents due to the difference of the diffu¬ 
sion velocity of M currents. The characteristic features of 
the multi-component fluid are the cross-correlation effects. 
For instance, Xab in Eq. dSTT i. in Eq. dSSl) . and so on 
arose from the correlation between microscopic currents as¬ 
sociated with the different conserved charges. Such coeffi¬ 
cients generate the correlation between the diffusions asso¬ 
ciated with different conserved charges. It is also notewor¬ 
thy that we have checked the positivity of the entropy pro¬ 
duction rate and the Onsager’s reciprocal theorem are natu¬ 
rally satisfied thanks to our adequate derivation of the multi- 
component hydrodynamics. Since we have faithfully solve the 
multi-component Boltzmann equation, these confirmation are 
also guaranteed by the Boltzmann theory. Furthermore, it is 
worth emphasizing that the form of the resultant relativistic 
hydrodynamic equation is of the same form as that derived in 
Ref. 11^ in which the positivity of the enetropy production 
rate and the reciprocal relations are imposed a priori in the 
derivation. It strongly suggests that the form of the reactive 
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multi-component relativistic equation has been now estab¬ 
lished, provided that the equation should satisfy the physically 
natural conditions, i.e., the positivity of the enetropy produc¬ 
tion rate and the reciprocal relations, irrespectively whether 
microscopic forms of the transport coefficients are given or 
not. 

Finally, we remark that the following studies are left as fu¬ 
ture works: We will take into account the particle creation and 
annihilation processes to analyze the realistic QGP hydrody¬ 
namic model. Such extension will make it possible to predict 
the observables quantitatively. 
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Appendix A: Detailed derivation of explicit form of excited 
modes 


In this section, we derive the expression of Eq. ( IA201 l. The 
initial value of the first-order equation belongs to the space 
which is spanned by as discussed in Ref. 11^ . 

Then, the calculation of L~^QqF^^'> can be reduced to that of 




with 


( 0 ) _ 1 


F^ ’ = 

k,Pk 


Pk ■ U 




M 


.4=1 


T 


(Al) 


(A2) 


Here, we have used Eq. (I38] l and the equilibrium one-body 
distribution function= [e^Pk-'^-'F.A’ik k'A)/T _ ^ 

To this end, we introduce the following quantities: 

^ I ipk • uY, f = 0, 1, • • • . (A3) 

Using these quantities, we have the metric 77 "^ = p^) as 

N 

n 




, k u ^ - Ml ^ 


k^l 


N 


■qP ^+3 = qMS P = ^ 


/C=l 


N 


n 


A+3 S+3 


= Y.qYq^MY 


(A4) 

(A5) 

(A 6 ) 




Then, we the inverse metric is written as 

ry-J = Au^^u'' + BA>^\ 
A-l-3 ~ '^A+3 fj. ~ CaU^ ^ 
'^A+3 B+3 ~ T^AB, 

where each coefficient satisfies 


N 


M N 


Y, Ml A - 1 + E E <ikMlCA = 0 , 

k^l A^lk^l 

N M N 

E qtMlA + E E <ikQkMlCB = 0 , 

k^l B^lk^l 

N M N 

E q^MlCB + E E <lhkMlVcB = 5ab. (A13) 


(A7) 

(A 8 ) 

(A9) 


(AlO) 

(All) 

(A12) 


fe=l C=l/c=l 

The inner products are evaluated as follows: 


N 


= E 


mlMl - Ml 


k=l L 

m’^Ml - Ml 


V T‘ 


1, 




M 


A^l 

N 




(A14) 


(^A+3^FW)=E'?fc 


aTtiIMI 


Ml -Ml 1 


k=l 


T 


V • u, (A15) 


Insetting the inverse metric (IA7l l- (IA91 l and the inner prod¬ 
ucts (IA141 l and (IA151 I into Eq. (lAlb . we have 

[QoF^X 


= H 


k,Pk 


M 


A=1 


P \-Lx7 

y / ' / ^ '^A,k,pk ( k- rp j + '^k,pk rp ! 


ly — 

h ^ T 


jP-V 

(A16) 


with (H, = (H, Jk-,Trk-'')k,p^/{pk-u) andllk,pk, 

JA,k,pk^ and are given by 

Ak,pk 




N 


( \ I n A — ^3 

+ (Pfe-w) Y^Aqk^^ -^ 

\ Z=1 


N M 

E E FABqYqt 

/=1 A,B=1 


AB — M 2 


--ml, (A17) 
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JA,k,pu = • «)) (A18) 

(A19) 

Finally, we obtain 'P as a linear combination of these bases 




[L-in] 


k,pk 


n 


+ 


(n,L-in) 

M 

A.B^l 

L _ -I fc,Pfc 

{7r'='^,L-i7rp£,) 


Jb,p 

(A20) 


with the initial value 

iioit = to', to) = Mo(fo) 

= C(fo)e'*°[/+ + C*(fo)e-“«t/-, (B7) 

where C'(fo) is a complex integration constant and C'*(fo) is 
its complex conjugate. U± = *(1, ±i) are eigenvectors of the 
matrix A satisfying AU± = ±iU±. 

The hrst-order equation with respect to e reads 

(^±-A^u = il-i:l)yol:iU+-U.). (B8) 

This equation is solved to be 

Mi(i;fo) = -^if+it',to)U+ + f-{t;to)U-), (B9) 


Here we have introduced the following nine vectors as mere 
coefficients of the basis vectors: 


n(CT;To), J^(ct;to), (a;To), (A21) 

which are identified as the bulk pressure, the heat flow, and 
the stress tensor, respectively. 


Appendix B: Solving the Van del Pol equation by RG method 

To show the prescription given in Sec. [Till we apply the RG 
method to solve the Van del Pol equation: 


with the dehnitions of the following quantities 

C*{to)-C{to)C*{to)^ 

2i ® 

-{t-to){Cito)-C{to)^C*{to))e" 


f-{t;to) = - 




4i 


3it 


2i 


2i 

-it-to)iC*{to)-C{to)C*{to)'^)e 


(BIO) 




(BID 


i + a; = e(l — (Bl) 

which is the oscillator with a limit cycle, which is a periodic 
behavior of the solution in the asymptotic regime. By rewrit¬ 
ing it in the form of simultaneous hrst-order differential equa¬ 
tion, we have 

(S 

where y = x and the following expressions are introduced 


An initial value reads: 

ui{t = to; to) = ui(to) 

= -^{f+{to;to)U+ + f.{to;to)U.). 

(B12) 

Here, we have added the zeroth-order solution so that the sec¬ 
ular terms vanish at f = to. One can hnd the detailed discus¬ 
sion in Ref. 

Summing up the solution up to 0(e), we have the perturba¬ 
tive solution and global solution 



First, we solve this equation perturbatively up to order e. We 
denote u{t) for the exact solution and u{t; to) for the pertur¬ 
bative solution, which is represented as a perturbation series: 

u{t;to) = Uo{t;to) + eui{t;to) + O(e^), (B4) 

with Ui = *{xi,yi). We take the initial value as the exact 
solution which is also expanded as 

u{to; to) = u{to) = uo(to) + eui(to) + O(e^). (B5) 

The zeroth-order equations with respect to e is solved to be 

(B6) 


u(t;to) = C(to)e'*U+ + C*(to)e-‘*U- 

-e^(f+(t;to)U+ + f.(t;to)U.), (B13) 

UG(t) = C'(De'*C/+ -b C*{t)e-^*U- 

-e^{f+{t;t)U+ + f.it;t)U.), (B14) 


respectively. 

Now, one can see that the RG equation: 


dtt 

dfo 


= 0 , 


(B15) 


to = t 


and the Van del Pol equation with the global solution inserted 


df 


- A] UG = € 


0 

(1 - xl;)yG 


uo{t;to) = C{to)e^*U+ + C*{to)e-^*U., 


(B16) 
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FIG. 2. The orbits of the numerical (exact) solution (a), the naive 
perturbative solution (b), and the RG-improved solution (c) of the 
Van del Pol equation (EB in the phase space [x,y) = (x,x). The 
initial value is set to (x,y) — (5.0, —1.4) commonly to the three 
solutions. 


indeed lead to the same differential equations up to 0(e): 

C(t) = l(C(t)-C(tfC*(t)). (B17) 

These equations have made the integration constant C(to) 
into the time-dependent variables. It is noteworthy that the 
right-hand side of the equation are 0(e), which indicates 
that we successfully identified the slow variable C(t) and ex¬ 
tracted the equations governing the slow dynamics, i.e., de¬ 
scribing how the radius of the orbit in the phase space gets 
relaxed to make a limit cycle. By parametrizing C(t) as 
C = (f?/2) exp(i0), the amplitude and angular motions may 
be extracted as 

= 0 = 0. (B18) 

Note that the former of Eq. (IBlSl l has a fixed point R = 2, 
i.e., the radius of the “limit cycle” toward which the orbit ap- 
praoches asymptotically: From FiglU we can see that while 
the permrbative solution (Fig.|2l(a)) behaves completely differ¬ 
ent from the numerical solution (Fig.|2l(b)) due to the secular- 
ity, the RG improved solution (Fig.|2l(c)) nicely reproduces the 
“limit cycle”, i.e., the characteristic property of the solution of 
the Van del Pol equation within the first-order perturbation. 


Appendix C: Detailed derivation of relaxation equation 

In this Appendix, we present a detailed derivation of the 
relaxation equation given by Eqs. (IC61l) - (IC631 l. The following 
calculation is based on a lengthy irreducible decomposition of 
various tensors, which are fully utilized in the derivation of 
Eqs. (IC21b - (IC28b and (IC44b - (IC60b . An explicit account of the 
tensor decompositions may be seen in a different document. 
We introduce the following “vectors”: 

(Cl) 

^i = {U, (C2) 

n/c.pfc J2b ^)bA I 

-rc’ T^/h ’ -2TTjj’ ^ ’ 

f 'y~'2 ^ ^ 

X, = UC9,—Y.^\AByp^,27ja^’'\, (C4) 

< = (C5) 

pk-u 

Then, Eq. ( l7Tb is converted into the following form: 

inn)-^[dr + ev ■ v][r'’(i + enL-^^xj)]) 
+ t^]^{L-^\B[L-^^Xo,L-^\k]) + 0(e"), (C6) 
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Expanding a term in the left-hand side we have 


e{L-^^\x^)X, 

+ e(L-^\ {nrr^ 


dr 


+ ev -V 


rnL-^<e)^j 


+ + 0(e")- (C7) 

The coefficients of the first term in the left-hand side and the 
first and second terms in the right-hand side of Eq. dCTll can 
be written as 


{L-^\X^) = 


-rr 2^ T’^jh ^ -2Tr] 

^ T'^/h ^ -2Tr) 


-TC 


^pcr 


\ -TC 2^ T’^/h ^ -2 Tt] 


1 0 0 

= I 0 hSAB^^^ 0 

0 0 


{L-^\L-^X^) = 




! {L ^n,L ^n) 11, 

-TC ^ T'^/h ^ -2Tr] 

- ^ — 


V 


-TC ^ Tyh 

{L-^TTP'',L-^fl) ^ {L-^ttP'',L-^J^) langleL-^T:P'^,L-^T^P^) 

-m 0 0 

0 -hrf^APP 0 
0 0 -t^AP'^p'^ 


-2Tr] 


(C8) 


(C9) 


^ Tyh ^ ^27^5 




-TC 

-TC 

~2fC 


c 




T‘^/h 


T^/h 


-2Tr] 

{L-^nP'',v°^L-^TTP'^) 


= I -7nA^^ 


-heijA-p 


—PB ^^ypoi 


0 


-i%AP^P^ 1 , 
0 


(CIO) 


r 


where we have introduced the relaxation times; 

Tn = ^(c-7c-'n)> 


(Cii) 




3T2 


M 

Y,{L-^J'X,L-^Jc,P){X-^)cb, (C12) 


c=i 


= 


lOTry 


{L-^ttP'^,L-^ttP''), 


(Cl 3) 


and the relaxation length: 


1 ^ 

^7 = E (A-')AB(C-'n, vPL-^Jb.p), (C14) 


B^l 




j/A _ 


lOTp 


{L ^J^,v''L ^TT^jy), 


(C15) 

(C16) 
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= — 


1 ^ 

^ (C17) 


5T2 




Then, the last term in the right-hand side of Eq. (IC71 t is written 


M 

= ^nnnll^ + ^ b^jjJ^JB.p + buTr-KT^^^T^pa, (C18) 

A,B = 1 




M 


— ^2 + b"jf^JB,pTTpp), 


(C19) 


B=1 


-i(L ^x^,L ^x''])V'jV’fe 


M 


= b^n.nn^^''+ J2 b^fj^^"^"JA,pJB,a 


A.B=1 


+ b^^^TTp^TTAo-, (C20) 


where the transport coefficients are defined by 

(L-^U, B[L-^Il, L-^U]) 


^nnn = — - 


2{Tcy 


,AB _ ^ {L-^Ii,B[L-^J^,L-^JD,p]) 

°njj = “ - 


C,D=1 

X (A“^)ca(A“^)b)b, 


6{T^/hy 


(C21) 


(C22) 


^ _ {L-^fl,B[L-^n^-,L-H^,]) 

bn.. - 10(2r77)2 ’ 

,AB (7.-iJ^,i?[L”ifl,L-iic.p]) 1 

bjnj = 2^ - otrr^\trr2/u\ -jCB, 


C=1 


OTOiT^/h) 


(C24) 

OJJtt = Z^ --(A JCB, 


C,B=1 


5(T2/fi)(2r,7) 


(C25) 

5(TC)(Tr?) ’ ^ ^ 

- Z^ 10(T2/fi)2 


C,B=1 

X (A“^)ca(A“^)b)b, 


(C27) 


_ (L-i7r^-,i7[L-i7r^,L-^7rA.]) 

""" “ (35/6)(2T7?)2 ■ ^ ^ 

Let us rewrite the second term in the right-hand side of 


— -f et; • V 

OT 


rrL-^x^Hj 


Eq. 0: 

+ ^ [rq/eqL-ix^])v>, ^ 

+ (l-^\ (/«qf 

+ e J2 {l-H\ {rr)-^vP^[rrL-^x^]y,^p^ 

(C29) 

The temporal derivative of T, ^jla/T, and are rewritten by 
using the balance equations up to the first order with respect 
to e, which correspond to the Euler equations: 


d 

—= -euAV ■ u + O(e^), 

d 

—e = —enhV ■ u + O(e^), 

OT 

AuM = eJ_V^P + 0(e2). 
OT nh 


(C30) 

(C31) 

(C32) 


These equations can be written as 


W T fi AI N „ 

/c=l S=lfe=l 


B k b) ^J■B 
Cio 


— — CTIaO H“ ); 

k^l B=1 k^l 

= -enhe + Oie"^), 

a 1 rp M N 

1 V^r + E E 


dT 


A=1 k=l 


T 


(C33) 

(C34) 

O(e^), 

(C35) 


where the definitions of a\, a\, and ^re given by Eq. (IA3]| . 
Let us define a matrix £ by writing Eqs. (IC331 l and (IC34I) as 


M 


^A0-^T+ ^ £ab-^^ - e6i + C>(e2), 

B—\ 


M 


b=i 


(C36) 

(C37) 


















16 


Therefore, Eqs. (IC33t - (IC331 l can be written as 


^T = ae + 0{e^), (C38) 

OT 

^^ = elAe + Oie^), (C39) 

= + V'^^ + 0(e2), (C40) 

A=1 I 


where we have defined the coefficients as I = X]s(^~^)ob 
and 2a = for notational simplicity. Then, 

Eq. ( IC29I ) takes the following form: 


L-ifi, inn)-^ 

M 


— +evS7 

OT 


nnL-^x^% 


= e 


M 


IVl IVl 

Knnll^ + 

^nj p~^ 

* - A,B=1 


per 


A^l 




Ji)^ 

^jn 


M M M M 

nv'^T + + Y + Y 

B=1 B=1 B=1 B=1 


M 

+ pT + Y 

B=1 


(nnr^ 


d 

— +evV 

OT 


nrL-Y^)xj 


= e 


(C41) 


(C42) 


M M 

^ KiYjA^T+ Y ^ ^ (C43) 

yt=l A.B=1 


where we have used V^Uv = a^v + Afj,^9/3 + oo^i, with the vorticity uj^v = iV^Uv — VvU^)l2, and the transport coefficients 
are defined as follows: 


'tnn = 


L-'n, {rrr 

AP!- 


\d d ^ a d 

A=1 J 


jeqjeq^-ljj 

~c 


(IM ^ (^eqjeq)-l 




n 


{2)AB _ 

nj — 

^IItt = 

.(1)^ = 
^jn — 

{2)AB _ 

,/n — 

Jl)^S _ 
'-JJ — 

MAb 
'‘JJ 


3 

AP 


d 1 a 


‘ dT T dui^ 


M 

E 


yeqjeq^-lj^^ 

-(A jsA 


P B=1 


T'^/h 


3 

^pvpa 


(l-^% innr^ 


d tiaT d 


nh 


-, M 

E 

C=1 


r2A//i ^ 


,-eq,-i, 9 nnL-H 


(^L-'ti.{f‘>n)-'v,— 


pa 


-2Tr] 




1 
3 


d Id 


/eq/eq^-in 


‘ dT T dui^ 


d hbT d 


-TC 
/eq/eq^-lfl 


a a 

i— + Yiv 


^ d^^ nh du^^ 
M 


-TC 


l^pvpa 


Qrp ^ QPD. - 

^ D=1 ^ T ] c=l 

M 


M 




T^X/h 


CB 




(C44) 

(C45) 

(C46) 

(C47) 

(C48) 

(C49) 

(C50) 

(C51) 
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_(3)AB_ d ^ /eq/eqL-lJc_ ^ 


l^fJ-vptJ 


Am ^ 

- — 


(L-^jA,^,{nn)-^ 


C=1 

d Id 

Vv-A^ + 




K 




Jtt 5 

^fii^pa- 

^ttII = Z 


(^L-^jA,^,{nn)-^ 


dT T du’' 


d ubT d 


d^ nh du'' 


-2Tr] 

nr^L-^p 

-2Tr] 


(l-1 V, inn) 


(l-Itt^., inn)-^ 


Jeq^-l 9 /®'1/®‘1L-In 


M 

E 

-■ S=1 


nnL-^jB,.,,.i, 

T^X/h 


^pvpcr 


„(2) = _ 

5 




= ±A>^mAP-\n^\s{L-^np,,{r^n)~"vc. 


"dm -TC 
d Id 

Vp-^ + y ^ 

A Z^A 

'^hdup 
_ 9 A T- 9 1 u 9 

A=i J J 

i 9 nnLnp 


/eq/eq^-V, 

T^X/h 


yeqjeq^-1^^ 


per 


-2T?7 


9u^ -2r7; 

9 nnLnp 

duP -2Tf] 


(C52) 

(C53) 

(C54) 

(C55) 

(C56) 

(C57) 

(C58) 

(C59) 

(C60) 


where fl^’^P'^ = (A^^’A'^'^ — A'^p)/2 is the projection operator onto an antisymmetric tensor of rank two. Substituting the 

above equations into Eq. (IC7]i, we arrive at the explicit form of the relaxation equations; 


9 


M 


en = -eC9 - ern^n + £2 I - ^ • Ja 


M 


M 


A^l 


A,B^1 


+ Knnn6» + A ^^A^A.pV^T + A + ^nTi-TTpo-o-' 


M \ 

+ ^Hnnll^ + A ^nJJ'^AB,p + pa j, 

A.B=1 / 


(C61) 


M 


M 


9 


eJ^ = eJ2 - £ A rf^^^^^JB,P + 

B=1 \ 


B=1 


M 


M 


M 


S = 1 


+ K^AnV'^T + A + A + A 

M 

^Vpr+AA^A'^^V p^ 

M \ 

+ A + A ^A'^B.pTT^'" ], 




M 


(C62) 




S=1 


M 


eTT 


/it/ 


= e2r,m'' - er^Af^'^P^—TTp^ - 


A 


V a—1 
M 


M M 

+ K^num'" + A 4AAv">r + A A^ 


A^l 

M 


AM^l 




(C63) 


A,B = 1 
































18 


After setting e = 1, we find that these equations become Eqs. (IC611 l- (IC631 l. 
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Note added: In the course of finalizing the present 
manuscript, we got aware of the paper by A. Jaiswal, B. 
Friman, and K. Redlich (arXiv; 1507.02849), in which the 
second-order hydrodynamics for the relativistic fluid com¬ 
posed of quarks, anti-quarks, and gluons is derived from the 
Boltzmann equation within the relaxation-time approxima¬ 
tion. It would be interesting to compare their results with ours 
to see how the simple relaxation-time approximation is useful 
and valid. 







